A role for the S4-domain containing protein YlmH in ribosome-associated quality control in Bacillus subtilis

Abstract Ribosomes trapped on mRNAs during protein synthesis need to be rescued for the cell to survive. The most ubiquitous bacterial ribosome rescue pathway is trans-translation mediated by tmRNA and SmpB. Genetic inactivation of trans-translation can be lethal, unless ribosomes are rescued by ArfA or ArfB alternative rescue factors or the ribosome-associated quality control (RQC) system, which in Bacillus subtilis involves MutS2, RqcH, RqcP and Pth. Using transposon sequencing in a trans-translation-incompetent B. subtilis strain we identify a poorly characterized S4-domain-containing protein YlmH as a novel potential RQC factor. Cryo-EM structures reveal that YlmH binds peptidyl-tRNA-50S complexes in a position analogous to that of S4-domain-containing protein RqcP, and that, similarly to RqcP, YlmH can co-habit with RqcH. Consistently, we show that YlmH can assume the role of RqcP in RQC by facilitating the addition of poly-alanine tails to truncated nascent polypeptides. While in B. subtilis the function of YlmH is redundant with RqcP, our taxonomic analysis reveals that in multiple bacterial phyla RqcP is absent, while YlmH and RqcH are present, suggesting that in these species YlmH plays a central role in the RQC.


Introduction
In all cells, protein synthesis on the ribosome is a multistep process that culminates in release of the finished protein (translation termination) and separation of the ribosome from the mRNA (ribosome recycling) ( 1 ,2 ).Ribosomal stalling on damaged or truncated mRNA molecules is harmful due to sequestration of ribosomes and production of potentially cytotoxic truncated protein products.Therefore, all domains of life have evolved diverse ribosome rescue pathways that recycle such stalled ribosomal complexes (3)(4)(5)(6)(7).The first bacterial ribosome rescue system to be discovered was the transtranslation system ( 8 ).The RNA component of the system, transfer-messenger RNA (tmRNA; encoded by the gene ssrA, 'small stable RNA A'), serves as a template encoding the Cterminal SsrA peptide tag (AANDENYALAA in Esc heric hia coli ) followed by a stop codon ( 9 ).Acting in concert with canonical bacterial translation termination machinery, the trans-translation system rescues ribosomes stalled on mRNAs lacking an in-frame stop codon, while simultaneously marking the incomplete polypeptides with a C-terminal degron tag for destruction by the ClpXP protease ( 10 ,11 ).The key element of the SsrA degron is the C-terminal Ala-Ala-COO-moiety ( 12 ,13 ) that is directly recognised by the AAA + unfoldase ClpX ( 10 ).
Despite the trans-translation system performing an essential function in bacterial cells, in most bacteria, including Bacillus subtilis and E. coli , ΔssrA strains tend to display only relatively minor phenotypes, such as increased susceptibility to antibiotics causing ribosomal miscoding (14)(15)(16) as well as heat sensitivity ( 15 ,17 ).This can be explained by the high degree of functional redundancy of trans-translation with additional ribosome rescue systems present in these organisms.For example, alternative rescue factors (Arfs) catalyse the release of the aberrant polypeptide from a stalled ribosomal complex, but do so without adding a degron tag (18)(19)(20).These systems either recognise the stalled ribosome and recruit canonical class-1 release factors RF1 and RF2 to cleave the stalled polypeptide (as exemplified by E. coli ArfA ( 18 ) and B. subtilis BrfA ( 19 )) or act directly as a peptidyl-tRNA hydrolase (as exemplified by E. coli ArfB ( 21 )).Functionally, both ssrA and arfA genes are dispensable when disrupted separately in E. coli , but simultaneous ΔssrA ΔarfA double-deletion results in synthetic lethality ( 22 ), which can be suppressed by overexpression of ArfB ( 21 ,23 ).Similarly, simultaneous loss of trans-translation and BrfA-mediated ribosomal rescue in B. subtilis results in a severe synthetic growth defect ( 19 ).These strong genetic interactions underscore the functional complementarity and redundancy between trans-translation and Arfmediated ribosomal rescue.
Trans-translation and Arf-mediated rescue systems do not operate on cytoplasmic eukaryotic ribosomes, making these bacterial systems attractive targets for novel antimicrobial agents ( 24 ,25 ).Instead, eukaryotes possess the so called ribosome-associated quality control (RQC) pathway that performs a similar function ( 7 ,26 ).Recently, a homologous RQC pathway has been identified in a diverse set of bacterial species, while being notably absent in E. coli ( 7 ,27 ).Similarly to Arf-encoding genes, bacterial RQC genes also genetically interact with ssrA .Simultaneous genetic disruption of both trans-translation and RQC systems in B. subtilis results in a growth defect at 37 • C (which is further exacerbated at elevated temperature), as well as increased sensitivity to translation-targeting antibiotics (28)(29)(30).The mechanism of RQC has been recently extensively characterized in B. subtilis (28)(29)(30)(31)(32)(33)(34).Like trans-translation, the RQC system also extends stalled polypeptides with C-terminal degron tags ( 32 ,35 ).
The RQC process can be divided into initiation (splitting of the stalled ribosome), elongation (extension of the nascent polypeptide with a degron tag) and termination (release of the tagged protein) phases.Both in eukaryotes and bacteria, ribosomal collisions serve as a substrate that is efficiently recognized by the RQC machinery ( 31 , 36 , 37 ).In B. subtilis , rescue of stalled and collided ribosomes (disomes) is initiated by the dimeric MutS2 / RqcU ATPase that splits the leading (stalled) ribosome into subunits, thereby generating a free 30S subunit and a 50S-nascent chain complex, 50S-NCC ( 31 ).Recruitment of B. subtilis MutS2 to stalled disomes is mediated by its KOW (Kyprides, Ouzounis, Woese) and Smr (small MutSrelated) domains ( 33 ).In E. coli , the SmrB nuclease cleaves the mRNA on the collided ribosomes ( 38 ), which is similar to the Cue2 and NONU-1 Smr RNases that perform analogous functions in yeast ( 39 ) and Caenorhabditis elegans ( 40 ), respectively.By contrast, the Smr domain of B. subtilis MutS2 seemingly lacks this RNase activity ( 33 ).
Regardless, the 50S-NCCs generated by MutS2 serve as a substrate for RQC elongation which in bacteria is catalysed by two factors: RqcH and RqcP ( 28-30 ,32 ).RqcH is homologous and functionally analogous to eukaryotic Rqc2 / NEMF (Nuclear Export Mediator Factor).Yeast Rqc2 catalyses the addition of C-terminal Ala and Thr tracts, resulting in socalled 'CAT tails' ( 41 ), whereas mammalian CAT tails were shown to contain Ala as well as Thr , Tyr , Gly, Asp and Asn ( 42 ).A critical step of the eukaryotic CAT tailing mechanism is the stabilization of P-site tRNA by eIF5A which enables the subsequent coordination of an A-site tRNA by Rqc2 / NEMF to promote peptide bond formation ( 43 ).Instead of the eIF5A homologue EF-P, bacterial RQC elongation involves RqcP, a member of the S4 RNA-binding domain protein family ( 44 ).Acting together with RqcH, RqcP catalyses the synthesis of poly-alanine tails ( 28 ,29 ).While RqcP is a close homologue of heat shock protein 15 (Hsp15), Hsp15 is an E. coli ribosome rescue factor that stabilizes P-tRNAs on stalled 50S subunits ( 45 ).However, it cannot functionally substitute RqcP in B. subtilis ( 28 ,30 ), despite the S4 domains of the two factors forming similar interactions with the P-site tRNA of the 50Snascent chain complexes (50S-NCC) ( 45 ).
The final step of bacterial RQC, release of the degrontagged protein, is the least well-understood stage of the process.In eukaryotes, a dedicated tRNA endonuclease Vms1 / ANKZF1 cleaves the CCA-3 -end of the P-site tRNA, thereby releasing the tagged polypeptide chain from the ribosome ( 46 ,47 ).By contrast, no dedicated RQC termination factor has been identified in bacteria, rather, a recent study by Svetlov et al. showed that the peptidyl-tRNA hydrolase (Pth) can release the poly-alanine-tagged nascent chains in B. subtilis ( 34 ).
By taking advantage of the strong genetic interactions between trans-translation and other ribosome rescue systems, we have employed transposon sequencing (Tn-Seq) to identify potential novel rescue factors in B. subtilis .In addition to identification of the well-characterised ArfA-type rescue factor BrfA ( 19 ), as well as the RQC elongation factors RqcP and RqcH ( 28-30 ,32 ), our Tn-Seq screen led to the identification of YlmH, a poorly characterized S4-domain-containing protein.Cryo-EM reconstructions of ex vivo purified YlmH reveal that YlmH interacts with 50S-NCCs.YlmH has a binding site analogous to RqcP and, similarly to RqcP, can co-habit 50S-NCC with RqcH, which is strongly indicative of YlmH being a RQC factor.Using in vivo reporter assays, we demonstrate that overexpression of YlmH in an RQC-incompetent ΔrqcP strain can restore the functionality of RQC-mediated poly-alanine tailing, indicating that YlmH is functionally redundant with RqcP.Finally, our phylogenic analysis indicates that some YlmH-containing bacteria do not encode RqcP, suggesting that in these phyla YlmH is likely to play a central role in RQC-mediated poly-alanine tailing.

Gene neighborhood analysis
The ylmH genomic context of selected taxa was analysed with FlaG2s ( 51 ) ( https:// github.com/GCA-VH-lab/FlaGs2 ).ylmH accessions and their respective assemblies were used as queries in FlaGs2 with default parameters and 15 flanking genes.

Phylogenetic distribution analysis
Preparation of HMMs profiles : RqcH sequences were collected by PSI-BLAST using B. subtilis RqcH (accession: WP_003232089.1)as the first query against ref-seq_select_prot database.Sequences with the Query Cover threshold exceeding 50% were collected after three iterations and then divided into each phylum.Sequence alignments were generated for each subgroup using MAFFT FFT-NS-2 v7.490 ( 48 ) and the HMMs profiles were built with HMMER v3.4 hmmbuild ( 49 ).In a similar way, HMMs profiles for YlmH were prepared by sequences collected through a two-iteration PSI-BLAST search using B. subtilis YlmH (WP_003232140.1) as a query.For Hsp15 and RqcP, representative Hsp15 and RqcP homologues from a previous study (see Table S2 in ( 28 )) were used as input for hmmbuild.HMMs profiles for SmpB (PF01668) and other S4-related domains (PF00163, PF00472, PF00579, PF00849, PF01728, PF03462, PF13275) were downloaded from InterPro ( 52 ).
HMM search and curation : Homologues of YlmH, RqcH, RqcP, Hsp15 and SmpB were identified from predicted proteomes of each GTDB representative genomes, release 214.1 ( 53 ), by HMMER v3.4 ( 49 ) hmmsearch with the HMMs profiles.The GTDB representative genomes were downloaded from NCBI using NCBI dataset command-line tool v0.3.SmpB homologues were identified by selecting sequences with a bit-score for PF01668 SmpB domain exceeding the gathering cutoff value.To identify YlmH homologues, sequences with at least one of bit-scores for HMMs profiles for YlmH subgroups exceeding 100 were selected and aligned with MAFFT FFT-NS-2.Then, those having the N-terminal domains were selected as YlmH.To identify RqcH homologues, sequences with at least one of bit-scores for the HMM profiles for RqcH subgroups exceeding 80 were selected.After the MAFFT FFT-NS-2 alignment, sequences harbouring conserved Asp-Arg or similar motifs that are crucial for RqcH function ( 32 ) were selected and categorized as RqcH.RqcP and Hsp15 homologues were identified according the following procedure.Sequences with a bit-score for HMMs profiles for RqcP or Hsp15 exceeding 40 were selected.Subsequently, those with a bit-score for HMMs profiles for the following domains equal to or greater than the gathering cutoff value were excluded: the Nterminal Ribosomal_S4 (PF00163) domain, which is found within uS4 proteins, the N-terminal tRNA-synt_1b (PF00579) domain, which is found in TyrS and TyrZ, the C-terminal PseudoU_synth_2 (PF00849) domain, which is found in RluB, YhcT, YlyB and YtzG, and the C-terminal FtsJ (PF01728) domain, which is found in YqxC.In addition, those with a bit-score for the S4_2 (PF13275) ≥45, and those classified as YlmH homologues were excluded.Sequences harbouring highly conserved Arg2 and Asp4 ( B. subtilis numbering) were selected and among them, those possessing the K14-R15-R16 motif and having a distance from Arg2 (or the corresponding Arg) to the C-terminus ≤110 residues were classified as RqcP; the remaining were classified as Hsp15 (see Figure 3 D in ( Ref. (28) ).
Phylogenetic tree : Phylogenetic tree downloaded from GTDB ( 53 ) were pruned to retain only the genomes used in the analysis by APE ( 54 ) and visualized and annotated by iTol v. 6 ( 55 ).

Tn-Seq sampling, library preparation and analysis
Preparation of transposon-inserted cell library : The in vitro transposition experiment was performed as described previously ( 56 ).Briefly, the himar1 transposon harbouring a kanamycin resistance marker gene was transposed from the plasmid pKIG645 by MarC9 transposase ( 56 ) into purified chromosomal DNA derived from either a wild-type 168 strain or its ssrA deletion mutant derivative (VHB257).The transposon-inserted DNA, in which gaps were filled, was used for the transformation of both the 168 strain and the ΔssrA mutant.Transformants were selected on LB agar plate supplemented with 3 μg / ml kanamycin following an overnight incubation at 37 • C. 2 × 10 5 and 1.5 × 10 5 cells were separately pooled for the wild-type and ΔssrA strains, respectively, and stored as transposon-inserted cell libraries with 25% glycerol at -80 • C.
DNA libr ary prepar ation and deep DNA sequencing : DNA libraries were prepared as per previously established protocol with minor modifications ( 57 ).Briefly, cell libraries (3 × 10 7 -2 × 10 8 cells) were grown in 50 ml of LB medium supplemented with 3 μg / ml kanamycin at 37 • C until OD 600 of 0.5.The cells were lysed with Nucleobond buffer set iii (MA CHEREY-NA GEL) and the chromosomal DNA was purified with Sera-Mag magnetic beads (Cytiva).The resulting DNA sample was then digested by MmeI (NEB) and re-purified on magnetic beads.DNA adapters were ligated to MmeI-digested DNA fragments by T4 DNA ligase (Takara), and the ligation products were purified on magnetic beads.Subsequently, PCR amplification was performed using primers ADPT-Tnseq-PCRPrimer (5 -AATGATA CGGCGA CC ACCGAGATCTAC ACTCTTTCC CTA CA CGACGCTCTTCCGA TCT-3 ) and P1-M6-GA T-MmeI (5 -C AAGC AGAAGACGGC AT ACGAGA TAGAC CGGGGACTT ATCA TCCAACCTGT-3 ), and the resulting products were resolved by polyacrylamide gel electrophoresis.DNA fragments with length of ≈ 148 bp were eluted from the crushed gels and purified by Nucleospin Gel and PCR Cleanup column.The DNA libraries were sequenced on the NextSeq 500 sequencing system (Illumina).
Tn-seq data analysis: Reads were processed as per the MA-GenTA protocol ( 58 ) using an in-house script developed with SeqKit ( 59 ), and processed reads were aligned to the reference B. subtilis 168 genomic sequence (Assembly: ASM904v1) using Bowtie version 1.3.1 ( 60 ).Fitness calculations were performed in accordance with the MAGenTA protocol using custom scripts implemented in R ( 61 ).The data are provided as Supplementary Table S2 and the code is provided as Supple-mentaryCode.zip.

Bacterial strains and plasmids
Strains (and information regarding their sensitivity to lincomycin), plasmids, oligonucleotides and synthetic DNA sequences used in this study are provided in Supplementary Table S3 .

Growth assays
B. subtilis 168 wild-type and deletion strains were pre-grown on LB plates overnight at 30 • C. Fresh individual colonies were used to inoculate liquid LB medium cultures to starting OD 600 of 0.01 at 37 • C. When the cultures reached mid-log phase (OD 600 of ≈0.4), they were diluted to OD 600 of 0.1.The resultant cultures were used to prepare 10-to 10 5 -fold serial dilutions which were then spotted onto LB agar plates with or without 1 mM IPTG.The plates were scored after 18 h incubation at either 37 • C or 49 • C.

Sucrose gradient fractionation
The experiments were performed as described previously ( 62 ).Briefly, B. subtilis strains were pre-grown on Luria Broth (LB) plates overnight at 30 • C. Fresh individual colonies were used to inoculate 200 ml LB cultures.The cultures were grown until OD 600 of 0.8 at 37 • C, the cells were collected by centrifugation and dissolved in 0.5 ml of HEPES:Polymix buffer [5 mM Mg(OAc)2] ( 62 ).Cells were lysed using a FastPrep homogenizer (MP Biomedicals) and the resultant lysates were clarified by centrifugation.Ten A 260 units of each extract were loaded onto 10-35% (w / v) sucrose density gradients prepared in HEPES:Polymix buffer [5 mM Mg(OAc) 2 ].Gradients were resolved at 36 000 rpm for 3 h at 4 • C. Both preparation and fractionation of gradients was done using a Biocomp Gradient Station (BioComp Instruments), and A 280 was used as a readout during the fractionation.

MS analysis of tagging sites on the reporter protein
The experiments were performed as described previously ( 33 ).Strains expressing the ApdA reporter were grown in 100 ml of LB with 20 μM of bortezomib until OD 600 of 0.5.Cells were harvested by centrifugation and frozen at -80 • C. Cells were then thawed in 2 × Celln Lytic B cell lysis reagent (Sigma) and 0.2 mg / ml lysozyme for 10 min gently rocking.The lysate was clarified by centrifugation at 20 000 × g for 30 min at 4 • C. 50 ul of Strep-tactin Sepharose beads (IBA) were added to the supernatant, and then the samples were incubated for 1 h at 4 • C. The beads were washed four times 5 min each at 4 • C with IP wash buffer (20 mM Tris pH 8, 100 mM NH 4 Cl, 0.4% Triton X-100, 0.1% NP-40).Protein was eluted with elution buffer (20 mM Tris pH 8, 100 mM NH 4 Cl, 5 mM Desthiobiotin) by shaking at 4 • C for 1 h.Then, 36 ul of elution was reduced with 100 mM DTT in 100 mM triethylammonium bicarbonate (TEAB) buffer at 58 • C for 55 min and then the pH was adjusted to 8. The samples were alkylated with 200 mM iodoacetamide in 100 mM TEAB buffer in the dark at room temperature for 15 min.Proteins were pelleted and resuspended in 100 mM TEAB and proteolyzed with 13.8 ng / μl of LysC (Wyco) at 37 • C overnight.Peptides were desalted on Oasis u-HLB plates (Waters), eluted with 60% acetonitrile (ACN) / 0.1% trifluoracetic acid (TFA), dried and reconstituted with 2% ACN / 0.1% formic acid.
Desalted peptides cleaved by LysC were analysed by LC-MS / MS.Then peptides were separated by reverse-phase chromatography (2-90% ACN / 0.1% formic acid gradient over 85 min at a rate of 300 nl / min) on a 75 um × 150 mm ReproSIL-Pur-120-c18-AQ column (Dr Albin Maisch, Germany) using the nano-EasyLC 1200 system (Thermo).Eluting peptides were sprayed into an Orbitrap-Lumos_ETD mass spectrometer through a 1-um emitter tip (New Objective) at 2.4 kV.Scans were acquired within 350-1800 Da m / z targeting the truncated reporter.Precursor ions were individually isolated and were fragmented (MS / MS) using an HCD activation collision energy of 30.Precursor (fragment) ions were analysed with a 200 automatic gain control (AGC), 118 ms maximum injection time (IT) at 60 000 resolution with 6 cycles.The MS / MS spectra were processed with Proteome Discover v2.5 (Thermo) and were analysed with Mascot v2.8.2 (Matrix Science) using RefSeq2021_204_Bacillus.S and a database with peptides from the NanoLuc-AdpA-BleR reporter protein.Peptide identifications from Mascot searches were processed within the Proteome Discoverer-Percolator to identify peptides with a confidence threshold of a 1% false discovery rate, as determined by an auto-concatenated decoy database search.Proteomics results are provided as Supplementary Table S4 .

Cryo-EM sample preparation
Cryo-EM sample preparation was performed as established in previous work ( 64 ) Briefly, B. subtilis 168 cells VHB715 ( ylmH trpC2 ; ylmH :: spcR ) expressing YlmH-FLAG from pHT01 p43 wRBS ylmH-GS5-FLAG wt tGyrA were grown at 37 • C in Lysogeny-Broth (LB) medium (Roth) supplemented with 5 μg / ml chloramphenicol and shaking at 145 rpm until OD 600 of 1.5.Cells were collected in 25 mM HEPES-KOH pH 7.5, 100 mM potassium acetate, 15 mM magnesium acetate, 0.1% (v / v) NP-40 and 0.5 mM Tris carboxy ethyl phosphene (TCEP) buffer supplemented with protease inhibitor cocktail (Roche), flash-frozen in liquid nitrogen and lysed under cryogenic conditions using a Retsch MM400 (Retsch).The lysate was cleared at 16 000 rpm in a JA-25.50 rotor (Beckman Coulter) for 15 min and incubated with anti-Flag M2 affinity beads (Merck) for 1.5 h at 4 • C on a turning wheel.After in-batch wash with 20 ml lysis buffer without protease inhibitors, the beads were transferred to a Mobicol column and washed with 4 ml of 25 mM HEPES-KOH, pH 7.5, 100 mM potassium acetate, 15 mM magnesium acetate, 0.01% DDM (dodecylmaltoside) and 0.5 mM TCEP buffer, after which the YlmH complexes were eluted using 0.2 mg ml −1 3 × Flag peptide (Sigma) in wash buffer.The complexes were then applied to grids for electron microscopy analysis or resolved on 4-12% NuP AGE SDS-P AGE gels (Invitrogen) by staining with Instant Blue (Expedeon).Sample volumes of 3.5 μl (8 OD 260 per ml) were applied to grids (Quantifoil, Cu, 300 mesh, R3 / 3 with 3 nm carbon) which had been freshly glow-discharged using a GloQube (Quorum Technologies) in negative charge mode at 25 mA for 90 s.Sample vitrification was performed using an ethane-propane mixture (37:63) in a Vitrobot Mark IV (Thermo Fisher Scientific), the chamber was set to 4 • C and 100% relative humidity and blotting was done for 3 s with no drain or wait time.

Cryo-EM data collection
Data were collected in an automated manner using EPU v.2.6 on a cold-FEG fringe-free Titan Krios G4 (Thermo Fisher Scientific) transmission electron microscope operating at 300 kV equipped with a Falcon 4 direct electron detector (Thermo Fisher Scientific).The camera was operated in electron counting mode and data were collected at a magnification of 96 000 × with the nominal pixel size of 0.83 Å and a nominal defocus range of −0.4 to −0.9 μm.A total of 25 683 micrographs in EER format were collected with 5.14 e / px / s and 5.31 s of exposure (corresponding to a total dose of 40 e / Å 2 on the specimen).No statistical methods were used to predetermine the sample size.The sample size was selected on the basis of a 3-day data collection, which was chosen to obtain a sufficient number of particles for data processing.

Cryo-EM data processing
Processing was performed using RELION 4.0.1 ( 65 ,66 ).The pixel size for processing was adjusted to 0.8 Å from the nominal 0.83 Å during data collection owing to best correlation with published ribosome models at this pixel size.Movie frames were aligned with MotionCor2 ( 67 ) using 5 × 5 patches followed by CTF estimation of the resulting micrographs using CTFFIND4 ( 68 ) using power spectra from the MotionCor run.The CTF fits were used to remove outlier micrographs with estimated resolutions greater than 15 Å, which retained 25457 micrographs.crYOLO 1.8.3 ( 69 ,70 ) with its general model (gmodel_phosnet_202005_N63_c17.h5) was used for particle picking, which resulted in 2 550 984 particles.These were extracted in a box size of 64 px at a pixel size of 4.8 Å and subjected to 2D classification.
After 2D classification, 2 391 525 particles resembling 50S subunits were selected and used for a first 3D auto-refinement to centre all particles for further refinement steps.An empty mature 50S subunit was used as reference, with the initial volume being generated from PDB ID 6HA8, ( Ref. (71) .Afterwards, particles were extracted with re-centring from the previous Refine3D-job at a box size of 128 px and a pixel size of 2.4 Å.The particles were aligned into a 3D volume using the output of the initial Refine3D-job as a reference (rescaled to the new box and pixel sizes).From these aligned particles, 3D classification was performed without further angular sampling.Particle sorting was performed according to Supplementary Figure S1 .Particles for final classes of the main YlmH-state and the YlmH-RqcH-state were re-extracted at a box size of 384 px with a pixel size of 0.8 Å and subjected to 3D auto-refinement.Particles for the main YlmH-state were further CTF-refined to correct for anisotropic magnification, trefoil and higher-order aberrations, defocus and astigmatism.Furthermore, particles for the main YlmH-state were subsequently subjected to Bayesian polishing followed by another round of CTF refinements.After these procedures, the final volumes were generated by 3D auto-refinement and postprocessing in RELION.

Molecular model building
The initial model for the 50S subunit of the main YlmHbound state was generated based on a published B. subtilis 70S structure (PDB ID: 6HA8) (71) .For YlmH, an AlphaFold ( 72 ,73 ) model was used and rigid-body fitted into the density using ChimeraX ( 74 ,75 ).Afterwards, the model was manually adjusted in Coot ( 76 ,77 ).Model refinement was performed using REFMAC5 as implemented in Servalcat ( 78 ) using pro-

B. subtilis S4-domain containing protein YlmH is a potential RQC factor
We have mounted a search for novel B. subtilis factors involved in ribosomal rescue by performing a Tn-Seq experiment with the ssrA (VHB257) strain as well as an isogenic parental 168 wild-type B. subtilis strain.The relative difference in the transposon insertion frequency between the two libraries reflects the difference in the insertion's effect on bacterial fitness in the two genetic backgrounds ( 83 ).The screen reliably identified genes encoding known ribosome rescue factors, including the release factor-dependent alternative ribosome rescue factor brfA ( 19 ) as well as the RQC factors rqcH and rqcP (Figure 1 A).Furthermore, a decreased frequency of insertion was observed for the ylmH ORF, indicative of YlmH also playing a role in ribosomal rescue.The experimentally unexplored ylmH gene is part of the ylm operon that is located downstream of the dcw ( division and cell wall ) ( 84 ) gene cluster ( 85 ).Importantly, just as RqcP and Hsp15, YlmH possesses a S4 RNA-binding domain ( 44 ), raising a possibility of a translation-related function.
We genetically probed the involvement of YlmH in ribosomal rescue by testing for a synthetic growth defect upon ylmH disruption with a spectinomycin resistant gene cassette ( spcR ) in the ssrA background.The growth phenotypes associated with ylmH and ylmH ssrA mutations were compared with that associated with rqcH and rqcP deletions.While the rqcH loss is not associated with a strong phenotype in B. subtilis under non-stress conditions, the ssrA rqcH double-deletion strain displays a growth defect at 37 • C and at elevated temperature (49 • C, mild heat shock conditions) the growth is abrogated completely ( 32 ) (Figure 1 B).The ssrA rqcP double-deletion results in a milder synthetic growth defect, with the effect being more pronounced at 49 • C. While the ylmH B. subtilis strain displays no growth defect at either 37 • C or 49 • C, the ssrA ylmH double-deletion phenocopies the ssrA rqcH strain, providing further support for the involvement of YlmH in RQC.
B. subtilis RQC factors are recruited to the 50S subunit upon challenge with translation-targeting antibiotics ( 29 , 30 , 32 ).Therefore, we were interested to assess whether YlmH also co-migrates with the 50S subunit under antibiotic stress.To do so, we used sucrose gradient centrifugation and immunoblotting to monitor the potential recruitment of C-terminally FLAG-tagged YlmH to ribosomal particles under different antibiotic stress conditions (Figure 1 C).In the absence of antibiotics, YlmH appeared to be predominantly located at the top of the gradient and only minimal association with ribosomal particles was observed.By contrast, in the presence of chloramphenicol YlmH appeared to migrate with the 50S subunit, as observed previously for RqcP and RqcH ( 28-30 ,32 ).Moreover, this association appeared to be specific for chloramphenicol stress, since negligible association with ribosomal particles was observed with other antibiotic stresses, such as the presence of erythromycin or lincomycin (Figure 1 C).It is noteworthy that chloramphenicol is a context-dependent elongation inhibitor that forms extensive interactions with an alanine residue at the -1 position of the nascent chain (86)(87)(88).As bacterial RQC is mediated by formation of C-terminal poly-alanine tails ( 29 , 30 , 32 ), specific stabilization of YlmH by chloramphenicol could be indicative of the factor's physical association with poly-alaninecontaining ribosomal complexes.

Cryo-EM structures of ex vivo YlmH-50S-NCC
Given the apparent association of YlmH with 50S particles (Figure 1 C), we employed affinity chromatography from a B. subtilis ylmH strain expressing YlmH with a C-terminal FLA G-tag (YlmH-FLA G) to isolate native YlmH-ribosome complexes for structural analysis.We have earlier employed this experimental strategy for other ribosome binding factors ( 45 ) including RqcH and RqcP ( 28 , 30 ).SDS-PA GE analysis of the pulldown samples revealed that 50S subunit ribosomal proteins co-purified with YlmH-FLAG (Figure 2 A), which was confirmed by mass spectrometry (see Materials and Methods).Together with the sucrose density gradient analysis (Figure 1 C), these findings indicate that YlmH indeed interacts with 50S subunits, analogous to other RQC factors such as RqcP and RqcH ( 28-30 ,32 ).The YlmH-50S complexes were then subjected to single particle cryo-EM analysis (see Materials and Methods).While the initial 3D classification indicated that the majority of ribosomal particles were 50S subunits, little density for additional ligands (tRNAs) or factors (YlmH) was observed ( Supplementary Figure S1 A-C).Therefore, 3D classification was performed using a mask encompassing the A-, P-and E-sites on the 50S subunit ( Supplementary Figure S1 D).This yielded two classes with an additional density adjacent to 23S rRNA helix 69 (H69) that we assigned to YlmH.In class 1 (154237 particles, 6.4%), YlmH was well-resolved and accompanied by additional density for a P-site tRNA, whereas in the major class 2 (1585755 particles, 66.3%) YlmH was less resolved and there was no apparent density for any tRNAs ( Supplementary Figure S1 E).The minor class 3 (61275 particles, 2.7%) contained 70S particles lacking any additional density for YlmH, so this class was not considered further.Given the well-defined density for YlmH and the presence of P-site tRNA, class 1 was subjected to further rounds of 3D classification, eventually yielding a well-defined YlmH-50S state ('State 1' with 66113 particles, 2.7% of the total) ( Supplementary Figure S1 F,G) containing YlmH and P-site tRNA (Figure 2 B).State 1 could be refined to generate a structure of the YlmH-50S complex with an average resolution of 2.0 Å (Table 1 and Supplementary Figure S3 ).
The excellent quality of the cryo-EM map density for State 1, particularly within the core of the 50S subunit ( Supplementary Figure S3 ), enabled the identification of nine distinct modifications within the 23S rRNA ( Supplementary Figure S4 ).This encompassed six modifications that have been visualized in the 23S rRNA of the E. coli 70S ribosome ( 89 ), namely, the four base methylations 5MU794 (Ec747), 5MU1968 (Ec1939), 2MG2474 (Ec2445), 2MA2532 (Ec2503); one ribose methylation at OMG2280 (Ec2251) and the dihydrouridine DHU2478 (Ec2449).In addition, we observe three modifications that have not been observed previously in E. coli or T. thermophilus 70S ribosomes ( 89-91 ), including two base modifications, namely, 5MU620 and 7MG2603, as well as one ribose methylation, OMG2582 ( Supplementary Figure S4 ).G2582 is present in the A-loop at the peptidyltransferase center and the 2 -O-methylation was shown recently to be mediated by the methyltransferase RlmP (formerly YsgA) ( 92 ).It is interesting that G2582 (Ec2553) neighbors U2581 (EcU2552), which is 2 -O-methylated in other organisms, such as E. coli , that lack the G2582 (Ec2553) modification, suggesting the presence, rather than the exact position, of a modification in the A-loop is important for function.

YlmH binds analogous to RqcP on the 50S subunit
The additional density in State 1 could be unambiguously assigned to YlmH (Figure 2 B,C).Unlike B. subtilis RqcP that is 86 aa and comprises only a single S4 domain, B. subtilis YlmH is 257 aa and comprises three domains (Figure 2 C).In addition to the C-terminal S4 domain (residues 171-257), YlmH contains an N-terminal domain (NTD, residues 1-74) as well as central domain (residues 87-166).The central domain of YlmH contains an RNA recognition motif (RRM) with high structural homology to RRMs found in DEADbox helicases, such E. coli DbpA ( 93 ,94 ) ( Supplementary Figure S5 A,B).A Dali search ( 95 ) revealed that the NTD exhibits the highest structural similarity to the R3H domain found in an unpublished structure (PDB ID 3GKU) of a putative RNA-binding protein from Clostridium symbiosum , which is homologous to the YidC1-associated Jag / KhpB protein (BSU_41030) in B. subtilis ( Supplementary Figure S5 C).While the NTD and C-terminal S4 domain were generally well-resolved (2.5-3.5 Å; Supplementary Figure S3 ), enabling sidechains to be modelled with high confidence for the regions interacting with the rRNA, the central RRM domain linking the N-and C-terminal domains was less well-resolved (3.5-5 Å; Supplementary Figure S3 ), enabling only rigid body fitting of an AlphaFold homology model ( 72 ).Nevertheless, the conformation and arrangement of the NTD and RRM domains observed for YlmH in State 1 appears very similar to that seen in the unpublished X-ray structure (PDB ID 2FPH) of Streptococcus pneumoniae YlmH lacking the S4 domain and determined in the absence of the ribosome ( Supplementary Figure S5 D).As expected, the C-terminal S4 domain of YlmH is located similarly on the 50S subunit to RqcP ( 28 ,29 ), where it interacts with both the P-site tRNA and H69 of the 23S rRNA (Figure 2 D, E).This is particularly evident for the two α-helices within the S4 domain of YlmH that are oriented to establish multiple interactions with H69, as reported previously for RqcP ( 28 ,29 ).By contrast, the five-stranded β-sheet in the S4 domain is shifted by 15-20 Å, when comparing the relative positions between YlmH and RqcP (Figure 2 F).As a result, YlmH establishes additional interactions with the P-site tRNA from residues located in the loop between β-strands β4 and β5.While this interaction is not possible for RqcP, we note that RqcP contains a short αhelix (residues 68-73; Highlighted in Figure 2 E), which is absent in the S4 domain of YlmH, that contacts the anticodon stem loop of the P-site tRNA ( 28 ,29 ).
In State 1, we observed additional density within the polypeptide exit tunnel, which is better resolved in the upper region of the tunnel and becomes disordered, probably due to flexibility, towards the tunnel exit (Figure 2 G).The density for the nascent chain at the PTC is consistent with four alanine residues extending from the CCA-end of the P-site tRNA (Figure 2 H).We also observed density for the antibiotic chloramphenicol bound within the A-site of the PTC (Figure 2 G,H).Although initially unexpected, chloramphenicol was used to select for the YlmH-FLAG expression plasmid and was therefore present in the growth medium upon lysis and could also be potentially present at low concentrations within the cell during growth (see Materials and Methods).Thus, chloramphenicol could bind to the YlmH-50S complexes during purification, if it was not already bound to the complexes within the cell.As mentioned above, chloramphenicol is a context-specific inhibitor that binds stably to the A-site of the PTC when alanine (and to a lesser extent serine or threonine) is present in the -1 position of the nascent chain ( 86 ), providing an explanation for the observed co-migration of YlmH with 50S subunits on sucrose gradients in the presence of this antibiotic (Figure 1 C).We also note that the interaction between the -1 alanine and chloramphenicol observed in the State 1 is near-identical to that observed in complexes programmed with alanine in the -1 position ( 88 ) ( Supplementary Figure S5 E, F).Thus, while we cannot exclude that the nascent chain represents a mixture of sequences (where the sidechains are averaged away and only the backbone density remains), we favor the interpretation that YlmH is bound to 50S subunits bearing poly-alanine tails, as expected for a factor involved in RQC ( 28-30 ,32 ).
The strong genetic and structural evidence for the involvement of YlmH in poly-alanine tailing during RQC, prompted us to further sub-sort our YlmH:P-tRNA:50S complexes in search of additional functional states.This led to the identification of a small population of 50S subunits (0.2%), which we termed State 2 ( Supplementary Figure S1 H,I), that contained in addition to YlmH and P-tRNA, also A-tRNA and RqcH (Figure 2 I).Despite the limited particle number (3770 particles), we were nevertheless able to refine State 2 to an average resolution of 3.4 Å (Table 1 and Supplementary Figure S6 ).While the local resolution of the core of the 50S subunit reached towards 3.0 Å, the ligands were significantly less well-resolved, presumably due to flexibility ( Supplementary Figure S6 ).Nevertheless, the resolution (5-10 Å) was sufficient to dock homology models for YlmH, RqcH as well as A-and P-tRNAs (Figure 2 J), revealing a functional state with remarkable similarity to that reported previously for an RqcP:RqcH:50S-NC complex, which also contained both Aand P-site tRNAs ( 29 ) ( Supplementary Figure S5 G,H).This further supports the notion that YlmH can participate in polyalanine tailing during RQC in B. subtilis .

Importance of both N-and CTD of YlmH interaction with the 50S subunit
Conserved Arg2 and Arg16 residues within the N-terminal region of the S4 domain of RqcP mediate essential interactions with 23S rRNA, to the extent that alanine substitutions of either residue abrogate the factor's functionality and phenocopy the rqcP strain ( 30 ).These two arginine residues are also conserved in YlmH suggesting an analogous functional importance (Figure 3 A).Indeed, YlmH forms extensive contacts with H69 of the 23S rRNA through both the YlmHspecific NTD region (amino acids 1-74) and the S4 domain (amino acids 171-257) (Figure 3 B), which it shares with RqcP.The conservation of the functional residues in the S4 domain and extensive contacts observed in cryo-EM reconstructions raised the question as to whether the rRNA contacts pro- vided by the additional YlmH-specific CTD are also essential for functionality.To address this question, we first tested the functional importance of Arg182 (corresponding to Arg2 of B. subtilis RqcP), Arg195 (which could be the functional equivalent of Arg16 of B. subtilis RqcP) and Lys197 (used as a control; does not form contacts with rRNA) of YlmH in phenotypic assays by assessing the heat sensitivity of B. subtilis expressing mutant YlmH variants in the ssrA background.As expected, alanine substitutions of R182 and R195 phenocopy the ylmH deletion, whereas the K197A substi-tution had no effect on YlmH functionality, since the ylmH K197A ssrA strain grew similarly to isogenic ssrA at both 37 • C and 49 • C (Figure 3 C).This suggests that the conserved Arg182 and R195 that mediate the interactions of YlmH with H69 of the 23S rRNA (Figure 3 B) are likely to be critical for the factor's functionality, whereas the non-conserved Lys197 is unimportant.
Next, we tested the essentiality of the N-terminal domain of YlmH (amino acid residues 1-166) which is absent in RqcP.To do so, we expressed a YlmH variant lacking the N-domain ( ylmH 2-167) in the ylmH ssrA background under the control of IPTG-inducible P hy −spank promoter ( 96 ) (Figure 3 D).As a positive control, we overexpressed the wildtype full-length YlmH, and as a negative control, we expressed the inactive YlmH R182A variant.As expected, the IPTGinduced expression of wild-type YlmH but not the ylmH 2-167 or ylmH R182A mutant variants suppressed the relatively mild growth defect of ylmH ssrA B. subtilis at 37 • C. Similarly, under the mild heat shock conditions (49 • C), expression of the mutant variants could not rescue the synthetic lethality either, suggesting that both ylmH R182A and ylmH 2-167 are non-functional.Collectively, these results suggest that the N-terminal domain of YlmH is required for the factors' functionality.
Finally, we used sucrose gradient centrifugation and immunoblotting experiments to probe the recruitment of the Cterminally FLAG-tagged YlmH 2-167 and YlmH R182A variants to 50S subunits upon a chloramphenicol challenge (Figure 3 E).In contrast to the wild-type YlmH-FLAG, which co-migrates with 50S subunits in the presence of chloramphenicol (Figure 1 C), neither the R182A-substituted YlmH variant nor the N-terminally truncated YlmH 2-167 variant were detected in the 50S fraction (Figure 3 E).This suggests that the inability of these variants to rescue the growth phenotype of the ssrA rqcH double deletion strain (Figure 3 D) most likely stems from the inability of these variants to stably interact with the 50S subunit.Thus, we conclude that in addition to the interaction between the S4 domain of YlmH and H69, the N-terminal domain also contributes to the stable binding of YlmH to the 50S subunit.This interpretation is consistent with the direct interactions observed between N-terminal residues, such as Asp3 and Gln6, and 23S rRNA nucleotides located within H68 (Figure 3 B).

YlmH overexpression can rescue the growth defect of the ssrA rqcH strain
To address whether YlmH is functionally interchangeable with RqcP, we initially investigated whether overexpression of YlmH can rescue the growth phenotype observed when the ssrA rqcP double-deletion strain is subjected to heat stress at 49 • C (Figure 4 A).We also performed the analogous experiment using the ssrA rqcH double-deletion strain, although we did not expect YlmH to able to functionally substitute for RqcH.Indeed, no restoration of the growth defect at 49 • C was observed upon P hy-spank -driven expression of YlmH in the ssrA rqcH strain (Figure 4 A).Conversely, expression of YlmH restored the growth defect of the ssrA rqcP double-deletion strain.This result is reminiscent of ArfB overexpression being able to suppress the synthetic lethality of a ssrA arfA double deletion in E. coli ( 21 ,23 ).Finally, to test the specificity of the effect of YlmH overexpression, we probed whether overexpression of other S4 domain-containing B. subtilis proteins can suppress the temperature sensitivity of the ssrA rqcP strain.The P hy-spank -driven expression of none of the tested genes ( rluB , rpsD , tyrS , tyrZ , yaaA , yhcT , yjbO , ylyB , yqxC and ytzG ) suppressed the heat-sensitivity phenotype of ssrA rqcP cells ( Supplementary Figure S7 ), suggesting that the rescue effect of YlmH overexpression is specific.

YlmH overexpression rescues poly-alanine tailing in ΔrqcP cells
Taken together, our results indicate that YlmH can functionally substitute for RqcP, and therefore YlmH may also be involved in the elongation phase of RQC.To test this hypothesis, we used the recently developed RQC reporter 33 ) expressed under the control of a constitutive P veg promoter ( 97 ).The reporter encodes an N-terminally twin-Streptagged NanoLuc luciferase conjoined via a 31-residue ApdA stalling motif from Amycolatopsis japonica ( 98 ) with the Cterminally FLAG-tagged bleomycin resistance protein (BleR).A ClpP inhibitor bortezomib ( 99 ) can be used to stabilise the poly-alanine-tailed RQC products ( 32 ).As we have shown earlier, the α-NanoLuc signal increases in the rqcH strain due to the absence of efficient tagging and targeting of the stalled product for ClpXP-mediated degradation ( 33 ).The reporter was inserted into the amyE locus coding for a nonessential α-amylase ( amyE ::P vegstrep 2 -NanoLuc -apdA -bleR -FLAG ) of either the wild-type 168 B. subtilis strain or B. subtilis strains lacking RqcH ( rqcH ), RqcP ( rqcP ) or YlmH ( ylmH ).The cultures were grown in the presence or absence of 20 μM bortezomib and the expression of the reporter protein was visualized on immunoblots with α-NanoLuc antibodies (Figure 4 B).To compare the effects across conditions and strains, NanoLuc levels were quantified, and the signal in the untreated wild-type strain was set as 1.
In the wild-type strain, the α-NanoLuc signal increased 4.4fold for the bortezomib-treated culture (Figure 4 B), reflecting the failure of ClpP-driven degradation of the poly-alaninetailed nascent polypeptide.Consistent with the lack of ribosome rescue and poly-alanine tailing in RQC-deficient rqcP and rqcH B. subtilis cells, in these strains the α-NanoLuc signal is high (3.1-3.6-foldincrease over the wild-type level) and bortezomib-insensitive (3.5-fold wild-type level), By contrast, the ylmH strain behaves similarly to wild-type.In the absence of bortezomib, the α-NanoLuc signal was 1.1 of the wild-type level, and it increased upon the bortezomib treatment to 4.1-fold over the untreated level.These results suggest that in the ylmH strain efficient poly-alanine tailing is mediated by canonical RqcH and RqcP factors.
As our genetic experiments suggested that YlmH can compensate for the loss RqcP but not RqcH (Figure 4 A), we next employed the NanoLuc reporter to test the effects of YlmH overexpression in wild-type as well as in the rqcP and rqcH strains, both in the presence and absence of bortezomib (Figure 4 C).The inactive R182A-substituted variant of YlmH was used as a negative control, and YlmH was Cterminally FLAG-tagged to monitor its inducible expression.For the wild-type B. subtilis strain, the expression of neither the wild-type YlmH, nor the R182A-substituted YlmH-FLAG variant, had any effect on the weak α-NanoLuc signal in the absence of bortezomib (Figure 4 C, left).Furthermore, expression of neither of the factors had any effect on the increase (3.4-3.5-fold) in the α-NanoLuc signal observed in the presence of bortezomib.This indicates that YlmH expression has no dominant negative effect on the canonical RQC pathway.In the rqcH strain, the elevated NanoLuc signal was insensitive to both bortezomib and YlmH expression, indicative of the absence of poly-alanine degron tailing (Figure 4 C, middle); as expected, YlmH cannot functionally compensate for loss of RqcH.In contrast, expression of YlmH in the rqcP strain in the absence of bortezomib led to a decrease (1.4-fold) in the α-NanoLuc signal to near-wild-type levels, and in the presence of bortezomib the signal increased 4.1-fold, (Figure 4 C, right).Importantly, this effect was specific, since expression of the inactive R182A-substituted YlmH-FLAG in the rqcP strain had no effect on the NanoLuc signal, which remained equally elevated (4.1-5.1-fold)regardless the absence or presence of bortezomib.Collectively, these results suggest that, when overexpressed, YlmH can substitute for RqcP in the poly-alanine tailing reaction.
Finally, to directly assess the effect of YlmH on poly-alanine tailing, we used mass spectrometry to quantify the levels of alanine-tailing on the ApdA reporter (Figure 4 C,D).We expressed the reporter in the wild-type strain as well as ΔrqcH and ΔrqcP knockout strains overexpressing YlmH, either just the wild-type (in the ΔrqcH background) or both wild-type and R182A-substituted YlmH variants (in the ΔrqcP background).The strains were grown in the presence of bortezomib to prevent the degradation of truncated or Ala-tagged proteins, the reporter was immunoprecipitated via its Nterminal Strep 2 -tag, the eluted proteins were digested with lysyl endopeptidase (LysC) and the resulting peptides were analysed with liquid chromatography coupled tandem mass spectrometry (LC-MS / MS).Ribosomes stall on the ApdA sequence at the RAPP motif with the first Pro codon in the P site and the second Pro codon in the A site ( 98 ), and our previous work has shown that poly-alanine tails are added after the first Pro residue, yielding RAP conjugated polyalanine tails ( 33 ).In the wild-type strain, we observed Alatailed peptides at this site, whereas in the ΔrqcH strain, there was no detectable poly-alanine tailing despite the overexpression of YlmH; this strict requirement for RqcH for Alatailing is consistent with previous studies ( 7 , 30 , 32 ).However, in cells lacking RqcP but overexpressing YlmH, the number of Ala-tailed peptides was similar to that observed in wildtype cells.This result directly shows that YlmH overexpression can compensate for RqcP's role in poly-alanine tailing in RQC.Finally, when the YlmH R182A mutant was overexpressed in ΔrqcP cells, very few Ala-tailed peptides were detected, further supporting the functional importance of R182 of YlmH.

Discussion
In this study, we demonstrate that YlmH can participate in the RQC pathway and is involved in poly-alanine tailing, where it can functionally substitute for its homologue RqcP.Until now, little was known about the molecular function of YlmH, with existing reports exclusively focused on the physiological effects of ylmH gene loss.A large-scale B. subtilis gene deletion library analysis showed that, with a possible exception of growth on glycerol-ammonium minimal media, the loss of YlmH is not associated with a significant decrease in fitness ( 100 ).In S. pneumoniae disruption of the ylmH gene results in early cell lysis at the end of the exponential phase ( 101 ), hinting at the factor's role in stress tolerance.A more prominent phenotype was observed in Synechocystis sp.PCC 6803.In this Cyanobacterium, disruption of the sll1252 gene encoding the YlmH orthologue leads to sensitivity to photoinhibition as well as to Ca 2+ and Cl − depletion, again pointing towards stress response defects ( 102 ).As mentioned, YlmH is located in the ylm operon, located downstream of the division and cell wall ( dcw ) cluster ( 85 ).While in some organisms, the ylm operon remains part of the dcw cluster, suggesting a potential role in cell division, it is notable that in many others, including B. subtilis , the ylm operon is separate from the dcw cluster ( 85 ) (Figure 5 A, Supplementary Table S1 ).
The first indication that YlmH is involved in a ribosome rescue pathway comes from our unbiased transposon mutagenesis screen identifying genes that produce a synthetic lethal phenotype when disrupted simultaneously with transtranslation (Figure 1 A).In addition to YlmH, we identified the alternative rescue factor BrfA, which plays a back-up function for the trans-translation system ( 19 ), as well as RqcH and RqcP which are directly involved in the RQC pathway ( 28-30 ,32 ).Although we did not identify peptidyl-tRNA hydrolase (Pth) that was recently reported to participate in the termination phase of RQC ( 34 ), this is not unexpected because spoVC , the gene encoding Pth is essential in B. subtilis ( 103 ) and most bacteria ( 104 ,105 ).We could recapitulate the genetic interaction between ylmH and ssrA by demonstrating a slow growth phenotype of the ΔylmH ΔssrA double deletion strain at 37 • C and genetic lethality at 49 • C (Figure 1 B).Moreover, we could show that YlmH can bind to peptidyl-tRNA-50S subunits (Figure 2 B) with a binding site that is similar to that observed previously for RqcP ( 28-30 ) (Figure 2 D-F).Additionally, we observed YlmH bound to peptidyl-tRNA-50S subunits containing RqcH (Figure 2 I, J), supporting a potential involvement of YlmH in RQC-mediated poly-alanine tailing.Although deletion of YlmH did not lead to a defect in RQC-mediated poly-alanine tailing, as we observed for RqcH or RqcP (Figure 4 B), we could show that overexpression of YlmH can rescue RQC mediated poly-alanine tailing when RqcP is absent (Figure 4 C,D).The rescue ability of YlmH does not appear to arise simply because it has an S4 domain like RqcP, but rather appears to be specific for YlmH since overexpression of other S4-containing proteins, including RspD, T yrS, T yrZ, YaaA, YhcT, YjbO, YlyB, Y qxC, YtzG and RluB, does not suppress the growth defect of the ΔrqcP ΔssrA double deletion strain ( Supplementary Figure S7 ).
Interestingly, the growth defect of the ΔylmH ΔssrA strain was more dramatic than that observed for the ΔrqcP ΔssrA and ΔrqcH ΔssrA double-deletion strains (Figure B), raising the question as to whether YlmH has additional key functions in B. subtilis .Indeed, YlmH contains additional N-terminal domains that are not present in RqcP (Figure 2 C), although we note that our findings indicate that the N-terminal domain of YlmH makes an important contribution to the 50S subunit binding activity (Figure 3 B,D).Nevertheless, strong support for a second function of YlmH comes from the cryo-EM analysis of the ex vivo YlmH pull-downs where we observe that the majority of the particles (66%) are bound to 50S particles lacking P-site tRNA ( Supplementary Figure S1 ).Sub-sorting of these particles revealed only states that lack density for H68 and present a shifted H69 which is an arrangement commonly observed during late steps of 50S biogenesis ( 106 ).This suggests, that YlmH may have an additional function in a late stage of ribosome biogenesis, however, this will require a further dedicated analysis.If YlmH does play a role in ribosome biogenesis, then this function does not appear to be essential since no growth defect was observed in the ΔylmH strain, although such a role may account for the stronger growth defect in the ΔylmH ΔssrA deletion strain compared to ΔrqcP ΔssrA and ΔrqcH ΔssrA deletion strains (Figure 1 B).
Although YlmH may only play a critical role in RQC in B. subtilis under specific conditions, an analysis of the distribution of YlmH across more bacterial taxonomies indicates that YlmH is likely to be a dedicated component of the RQC machinery in some bacterial lineages (Figure 5 B,C and Supplementary Table S1 ).Indeed, a YlmH-RqcH co-distribution to the exclusion of RqcP is characteristic of Deinococcota (95% out of 133 genomes analysed), suggesting that in these bacteria the YlmH-RqcH pair alone is responsible for catalysing poly-alanine tailing.Similarly, half of the analysed cyanobacterial genomes, including the experimentally studied Synechocystis sp.PCC 6803, encode YlmH and RqcH, but no RqcP or Hsp15, indicative of RQC being driven by the YlmH-RqcH pair in these species.By contrast, Bacillota (formerly Firmicutes) genomes, typically encode YlmH, RqcP and RqcH (64% of 10426 genomes analysed) suggesting that the functional redundancy in RQC observed in B. subtilis may be typical for this phylum.Similarly, 85% out of 86 Fusobacteriota genomes encode YlmH, RqcP and RqcH.The fact that both homologues have been maintained throughout evolution of different lineages of Bacillota and Fusobacteriota indicates an important partitioning of functions.S4domain factors Hsp15 and RqcP are largely mutually exclusive in distribution, with RqcP strongly co-distributing with RqcH ( 28 ).Thermotogota, Synergistota, Zixibacteria, Cloacimonadota and Campylobacteria genomes predominantly encode the RqcH-RqcP pair but no YlmH (79%, 84%, 67%, 77% and 92%, respectively).While the sll1252 -disrupted Synechocystis sp.PCC 6803 strain that was characterized by Inoue-Kashino et al. ( 102 ) is likely to be RQC-deficient, this hypothesis awaits experimental validation.Finally, significant fractions of Armatimonadota, Spirochaetota, Myxococcota and Nitrospirota genomes (69%, 65%, 33% and 43%, respectively) encode RqcH but no RqcP or YlmH, suggesting that if these organisms possess a functional RQC pathway, some yet-to-be-identified factor performs the function of RqcP / YlmH.

Figure 1 .
Figure 1.B. subtilis YlmH is implicated in Ribosome Quality Control.( A ) Deletion of ylmH in the ssrA background is associated with a fitness loss.Transposon libraries were generated in 168 wild-type and ssrA (VHB257) B. subtilis strains.Sites of transposon insertion were identified by Illumina NextSeq sequencing and mapped onto the B. subtilis 168 reference genome.Independent insertions in the genomic region encompassing ssrA , brfA , rqcH , rqcP and ylmH genes as well as insertion's fitness effect are shown.( B ) Deletions of ylmH , rqcH and rqcP in the ssrA background cause a synthetic growth defect.10-fold serial dilutions were spotted onto LB agar plates and incubated for 18 h at either 37 • C (left; optimal temperature) or 49 • C (right; mild heat shock).( C ) Chloramphenicol treatment promotes YlmH recruitment to 50S.C-terminally FLAG-tagged YlmH (YlmH-FLAG; C-terminal FLAG tag connected via GS5 linker) was ectopically expressed in the ylmH background under control of the P h y -spank promotor (BCHT931).YlmH e xpression w as induced b y 30 μM IPTG.Follo wing a 20-min antibiotic treatment with either chloramphenicol (Cam, 5 μg / ml), erythrom y cin (Ery, 1 μg / ml) or lincom y cin (Lnc, 80 μg / ml), cellular lysates were resolved on sucrose gradients and fractions probed with anti-FLAG antibodies.The polysome profile shown at the top represents the untreated condition.The full set of polysome profiles for all antibiotic treatments are shown in Supplementary Figure S2 A.

Figure 2 .
Figure 2. Cryo-EM str uct ure of YlmH-50S complexes.( A ) SDS-PAGE analysis of control (-) and YlmH-FLAG coimmunoprecipitations revealing the presence of ribosomal proteins of the large subunit.The background band visible in both samples corresponds to enolase.( B ) Cryo-EM reconstruction of State 1 consisting of YlmH (orange), P-tRNA (blue) bound to the large 50S subunit (y ello w) with landmarks L1-stalk and central protuberance (CP) highlighted.( C ) Schematic of the domain str uct ure of YlmH with N-(purple), central RRM (red) and S4 (orange) domains indicated.L o w er panel shows the isolated cryo-EM density (transparent grey) for YlmH in the same orientation as in panel (B), together with the molecular model (coloured by domain).(D-F) comparison of the relative binding position of ( D) YlmH (orange), ( E ) RqcP (green, PDB ID 7AS8) ( 28 ) to P-tRNA (blue) and ( F ) superimposition of (D) and (E) with arrows indicating the shift in position of the β-sheet of YlmH compared to RqcP.The extra α-helix of RqcP is highlighted in (E).Alignments were made on the basis of the 23S rRNA.( G ) Transverse section and zoom of the polypeptide exit tunnel of State 1 re v ealing the presence of a polypeptide chain (red) attached to the P-site tRNA (blue) and chloramphenicol (y ello w) bound in the A-site of the 50S subunit (light blue).( H ) Isolated cryo-EM density (transparent grey) of the CCA-end of the P-tRNA (blue), chloramphenicol (yellow) and nascent chain (red, with poly-alanine model) with selected 23S rRNA nucleotides (grey).( I ) Cryo-EM reconstruction of State 2 consisting of YlmH (orange), P-tRNA (blue), A-tRNA (green) and RqcH (mustard) bound to the large 50S subunit (y ello w) with landmarks L1-stalk and central protuberance (CP) highlighted.( J ) Isolated cryo-EM density (transparent grey) and fitted molecular models for YlmH (orange), P-tRNA (blue), A-tRNA (green) and RqcH (mustard) from State 1.

3 .
Contacts with H69 23S rRNA and P-site tRNA are essential for YlmH functionality.( A ) Sequence alignment of YlmH homologues from representative RqcP-encoding and RqcP-lacking bacterial species.( B ) Interaction network between 23S rRNA H69 and YlmH.R182 and R195 are predicted to interact with H69, whereas K197 is not.( C ) Disruption of YlmH contacts with H69 abrogates YlmH functionality in vivo .Synthetic growth defects ylmH mutant variants in the ssrA background of were compared with ylmH deletion.10-fold serial dilutions were spotted onto LB agar plates and incubated for 18 h at 37 • C (left) or 49 • C (right). ( D ) R182A substitution and N-terminal truncation ( 2-187) of ylmH abrogate its functionality in in vivo.W ild-t ype and mut ant ylmH v ariants w ere e xpressed in ylmH ssrA back ground under the control of P h y -spank promotor . 1 0-fold serial dilutions were spotted onto LB agar plates with or without 1 mM IPTG, and incubated for 18 h at either 37 • C or 49 • C. ( E ) R182A and 2-187 YlmH variants compromised in H69 binding are not recruited to the 50S upon chloramphenicol challenge.Sucrose gradient sedimentation and anti-FLAG immunoblotting of YlmH R182A -FLAG (BCHT1327) and YlmH 2-187 -FLAG (BCHT932) expressed in the ylmH background, with or without 5 μg / ml chloramphenicol (Cam) treatment.The polysome profile shown at the top represents the untreated condition of the YlmH R182A variant.The full set of polysome profiles for all analysed gradients are shown in Supplementary Figure S2 B.

Figure 4 .
Figure 4. YlmH can substitute RqcP, but not RqcH, during poly-alanine t ailing .( A ) Ectopic expression of ylmH suppresses the growth defect of the rqcP ssrA but not of the rqcH ssrA strain.Strains with or without IPTG-inducible construct in either rqcH ssrA (VHB788, VHB787) or rqcP ssrA (VHB789, VHB600) background were serially diluted, spotted onto LB agar plates without 1 mM IPTG and incubated for 18 h at either 37 • C or 49 • C. ( B ) rqcH and rqcP but not ylmH strains are defective in RQC-mediated poly-alanine t agging .A 31-residue ApdA stalling motif from A. japonica was inserted between N-terminally twin-strep tagged NanoLuc and C-terminally FLAG tagged bleom y cin resistance protein (BleR).The reporter is expressed under the control of constitutive P veg prompter in wild-type (BCHT1294), rqcH (BCHT1295), rqcP (BCHT1296) and ylmH (BCHT1297) strains after treatment with 20 μM ClpP inhibitor bortezomib and detected with anti-NanoLuc antibodies.The SigA protein was used as a loading control.For quantification the anti-NanoLuc signal was normalized to the internal anti-SigA control, and results of three independent biological replicates are shown as mean ± SD. ( C ) YlmH but not inactive YlmH R182A can substitute RqcP in poly-alanine tailing reaction.Wild-type ylmH or ylmH R182A were ectopically expressed in wild-type (BCHT1133, BCHT1134; left ), rqcH (BCHT1135, BCHT1136; middle ) and rqcP (BCHT1137, BCHT1138; right ) strains, and accumulation of the NanoLuc reporter was assessed with or without 20 μM bortezomib treatment.Data quantification was done as for panel (B).The arrow indicates the position of C-terminally FLAG-tagged wildtype YlmH and YlmH R182A .( D ) Ov ere xpression of wild-type but not R182A YlmH rescues poly-alanine tailing in rqcP but not the rqcH strain.The Ala-tail at the stall site of the reporter in wild-type, rqcH , and rqcP strains expressing either wild-type ylmH or ylmH R182A was detected by LC-MS / MS.The peptide counts are normalized to a different peptide in the reporter protein.

Table 1 .
Cryo-EM data collection, refinement and validation statistics